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The pressure and temperature dependence of the phonon dispersion of 2H-NbSe2 is measured by 
inelastic X-ray scattering. A strong temperature dependent soft phonon mode, reminiscent of the 
Charge Density Wave (CDW), is found to persist up to a pressure as high as 16 GPa, far above the 
critical pressure at which the CDW disappears at 0 K. By using ab initio calculations beyond the 
harmonic approximation, we obtain an accurate, quantitative, description of the (P,T) dependence 
of the phonon spectrum. Our results show that the rapid destruction of the CDW under pressure 
is related to the zero mode vibrations - or quantum fluctuations - of the lattice renormalized by the 
anharmonic part of the lattice potential. The calculations also show that the low-energy longitudinal 
acoustic mode that drives the CDW transition barely contributes to superconductivity, explaining 
the insensitivity of the superconducting critical temperature to the CDW transition. 


The interplay between charge density wave order, i.e. 
a static modulation of the electronic density close to the 
Fermi level, and superconductivity has attracted much 
attention^’^. This is the central issue of a long standing 
debate in simple transition metal dichalcogenides with¬ 
out strong electronic correlations, such as 2 H-NbSe 2 ^’"’’. 
At Tc£uv= 33.5 K, 2H-NbSe2 undergoes a second order 
phase transition towards an incommensurate CDW phase 
which coexists with superconductivity below Tc = 7.2 K. 
The anisotropy and temperature dependence of the elec¬ 
tronic band structure has been widely studied by ARPES 
and STM measurements^^^^. The CDW is coupled to 
a periodic lattice distortion through a strong electron- 
phonon coupling. The transition is associated with a soft¬ 
ening of a longitudinal acoustic phonon mode as temper¬ 
ature is lowered to Tcdw^^- It was already noticed that 
the high order phonon fluctuations and strong electron- 
phonon interactions explain some of the key features of 
the formation of the CDW in this system^^d4 
specific case of NbSe 2 at ambient pressure, the reduction 
of the phonon lifetime through the electron-phonon cou¬ 
pling is a key ingredient on the temperature dependence 
of the phonon dispersion^®. 

At ambient pressure, published ab initio calculations 
reproduce successfully the lattice instability on a wide 
part of the Brillouin zone around the experimentally ob¬ 
served qcDW However, these calculations are 

carried out at the harmonic level, and therefore at T=0K. 
The temperature dependence can be qualitatively repro¬ 
duced using a Gaussian smearing of the Fermi surface 
which reduces the contribution of the electron-phonon 


interaction. Quantitatively however, the temperatures 
that would effectively yield this smearing are unphysi¬ 
cal, as they exceed the experimentally observed one by 
several orders of magnitude. 

The quantitative failure of the harmonic ab-initio cal¬ 
culations is even worse under pressure. Experimentally, 
by applying an hydrostatic pressure, Tcdw is driven to 
0 K and the CDW instability disappears above the crit¬ 
ical pressure Pcdw^ = 3.5 K)=4.6 GPa^®. Around 
this quantum phase transition superconductivity remains 
essentially unaffected, as shown on the (T,P) phase dia¬ 
gram in Fig. li9"20 However, harmonic ab initio calcula¬ 
tions fail to describe the ground state of NbSe 2 as a func¬ 
tion of pressure (see SI and Indeed, they predict (see 
Fig. 3 and SI) a CDW instability up to Pc = 14 GPa, a 
pressure that largely exceeds the experimentally observed 
one. This persistence of the CDW instability in the har¬ 
monic calculations is reminiscent of the effect of iso-valent 
and iso-electronic substitution of Se with lighter S. The 
harmonic calculation also predicts a CDW instability for 
NbS 2 that is not observed experimentally^^. 

Taken together, these facts indicate that for an ac¬ 
curate description of the ground state (even at T=0K) 
of these systems, calculations beyond the standard har¬ 
monic approximation must be carried out. The impact of 
anharmonicity has recently been highlighted in different 
systems. To mention only a few, the thermodynamical 
properties of the metal-insulator transition in strongly 
correlated V 02 ^®, the transport properties of thermoelec¬ 
tric PbTe^^ or the electron-phonon interaction above the 
Verwey transition in magnetite Fe 304 ^® are all strongly 
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FIG. 1. Phase diagram of NbSe 2 . Squares indicate the 
limit of the CDW phase measured in Ref.'^® with the 
addition of the point at 4.6 GPa from Ref.^®, triangles 
indicate the superconducting Tc measured in®®, and 
circles indicate the pressure and temperature of our 
IXS measurements. Inset : crystallographic structure 
of 2H-NbSe2®. 


influenced by anharmonicity. Even at low temperatures, 
where its effects are expected to be weaker, it must be 
taken into account to explain the pressure dependence of 
the crystal structure of Calcium®®, the inverse isotope ef¬ 
fect in superconducting PdH®^, or the phonon spectra at 
the proximity of the liquid-solid transition of Helium at 
OK®®. Finally, anharmonicity is crucial for understand¬ 
ing the ground state of a system when a quantum elec¬ 
tronic phase transition coupled to the lattice occurs, as 
in e.g. ferroelectrics®® or CDW compounds^®’^"^’®®’®®. In 
most cases however, these effects are computationally too 
expensive to be evaluated in standard first-principle cal¬ 
culations, and therefore often disregarded. 

Here, we report the investigation of the temperature 
and pressure dependence of the phonon spectra of NbSe 2 
by using high resolution inelastic X-ray scattering (IXS) 
from a crystal in a diamond anvil cell, up to pressures 
as large as 16 GPa at T = 20 K. We observe that the 
strong temperature dependence of the soft phonon mode 
is reduced by pressure, but still present even at 16 GPa. 

We show that the temperature and pressure depen¬ 
dence of the phonon spectra can be accurately accounted 
for when the effects of anharmonicity are explicitly in¬ 
cluded in the calculation. Anharmonic effects were 
calculated within the newly developed stochastic self- 
consistent harmonic approximation (SSCHA), yielding 
an unprecedented agreement with the experimentally de¬ 
termined temperature dependence of the phonon spectra 
at various pressures. We show that at low temperatures, 
for pressures between Pcdw and Pc, the CDW is de¬ 


stroyed by quantum fluctuations, which easily overcome 
the double-well lattice potential and, thus, are strongly 
affected by the anharmonic part of the potential. Surpris¬ 
ingly, we observe that anharmonic effects dominate the 
low-energy phonon spectrum of 2H-NbSe2 in a large pres¬ 
sure range, extending at least up to 16 GPa >> Pcdw 
and down to the lowest temperatures. Finally, we demon¬ 
strate that a large electron-phonon interaction, mostly 
due to optical modes rather than to the soft longitudinal 
acoustic mode which drives the CDW instability, con¬ 
tributes substantially to superconductivity in 2H-NbSe2, 
and naturally explains its insensitivity to the occurrence 
of a CDW. 

Single crystals of 2H-NbSe2 were grown using the va¬ 
por growth technique in sealed quartz tubes with iodine 
as a transport agent®. Typical dimensions of the crystals 
are 100 x 100 x 50 ^m® (a x 6 x c). Pressure was generated 
using a diamond anvil cell, loaded with helium as a pres¬ 
sure transmitting medium to ensure highly hydrostatic 
conditions. Temperature was lowered using a custom- 
designed ^He cryostat, and pressure was varied in-situ at 
low temperature using a helium-pressurized membrane. 
Each cell contained two rubies to monitor the pressure in- 
situ using the fluorescence technique. Measurements of 
the phonon dispersion of 2H-NbSe2 were taken between 
room temperature and 20 K, and for pressures ranging 
from 0.2 to 16 GPa. 

IXS measurements were carried out on beamline ID28 
at the ESRF. The x-ray beam was aligned along the c- 
axis of the crystal. We used the (9,9,9) reflection on the 
backscattering monochromator with an incident energy 
of 17.794 keV, and a corresponding energy resolution of 
1.3 meV HWHM (lorentzian fit of the elastic peak). The 
incident beam was focused by a multilayer mirror into a 
spot of 100x60/im (width x height). We used 20x60mm 
(width X height) slits at a 7 m distance giving a momen¬ 
tum resolution of 0.014 A"® (a* x c* « 2.1 x 0.50 A"®) in 
the (H,0,L) plane, and 0.042 A“® (6*xsin(6CP) « I.SA”®) 
in the perpendicular direction. We measured the phonon 
dispersion along the TM direction close to r 2 oo (*■ e- 
(2-h, 0, 0)), where both longitudinal optical and acoustic 
like soft phonon branches can be observed. According 
to the balance of the dynamical structure factor of the 2 
phonon modes and so of the inelastic scattering function 
S(q,a;), the observation of one of the two branches can 
be selected. For each transfer wave vector investigated, 
an energy scan up to 15meV was performed. 

All the spectra were fitted using the standard 
Levenberg-Marquadt algorithm, with the following free 
parameters : position and amplitude of the elastic peak; 
position, linewidth and amplitude of the phonon; and 
a constant background. We make the standard assump¬ 
tions®®’®® that the elastic peak is Lorentzian, and that the 
phonon lineshape can be modeled by a damped harmonic 
oscillator convoluted by a normalized Lorentzian with a 
width corresponding to the resolution of the detector. 

Ah initio calculations were performed within den¬ 
sity functional perturbation theory and the general- 
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ized gradient approximation^^ making use of the Quan¬ 
tum ESPRESSO^^ code. An ultrasoft pseudopoten¬ 
tial (norm-conserving) for Nb (Se), a 35 Ry energy cut¬ 
off, and a 24x24x8 mesh for the electronic integrations 
were used. An Hermitian-Gaussian smearing of 0.01 Ry 
was used. Harmonic dynamical matrices were calculated 
within linear response in a grid of 6x6x4 q points. The 
small smearing parameter and the dense electronic mesh 
used are needed to converge the phonon frequency of the 
longitudinal acoustic mode at the CDW wave vector (cf. 
SI). The stochastic self-consistent harmonic approxima¬ 
tion (SSCHA)^’^’^"^ calculations were performed using a 
3x3x1 supercell, yielding anharmonic dynamical matri¬ 
ces in the commensurate 3x3x1 grid. The difference be¬ 
tween the SSCHA dynamical matrices and the harmonic 
dynamical matrices in the 3x3x1 grid was Fourier inter¬ 
polated to the points of the finer 6x6x4 grid. Adding 
the harmonic dynamical matrices in the 6x6x4 grid to 
the result of the interpolation, the SSCHA dynamical 
matrices in the finer 6x6x4 mesh were obtained. The 
experimental lattice parameters^® with relaxed internal 
coordinates were used because in these conditions the 
results weakly depend on the exchange-correlation po¬ 
tential chosen. 

The electron-phonon calculations were performed by 
using maximally localized Wannier functions®®’®® and in¬ 
terpolation of the electron-phonon matrix elements as in 
Ref.®^. We used 14 Wannier functions (d-states of Nb and 
Se p 2 states) and a 6 x 6 x 4 electron-momentum grid. 
We used the Fourier interpolated SSCHA dynamical ma¬ 
trices in the electron-phonon coupling calculations. The 
average electron-phonon coupling was interpolated using 
electron and phonon momentum grids of 24 x 24 x 8 ran¬ 
domly displaced from the origin. 

The temperature dependence of the low energy phonon 
dispersions recorded at pressures of 3.5 and 16 GPa are 
shown in Fig. 2 (some of the corresponding raw data 
is shown in the supplementary materials). At T =20 
K and P = 3.5 GPa, the system is very close to the 
quantum critical regime reported in ref. 18. No particu¬ 
lar enhancement of the elastic line intensity close to the 
CDW wavevector qcDW ~ (1.67,0,0) is observed. On 
the other hand, a colossal temperature dependent soft¬ 
ening of the longitudinal acoustic phonon mode is seen : 
from ~ 9 to ^ 1.5 meV between room temperature and 20 
K. While this mode is expected to condense around ~ 12 
K®®, this temperature could not be reached at this pres¬ 
sure with the current experimental set-up. As we increase 
pressure to 16 GPa, the amplitude of the phonon soften¬ 
ing is significantly reduced, but remains sizeable as the 
mode’s energy almost decreases by a factor of 2 from 300 
to 20 K. So far, most of the theoretical and experimental 
work on NbSe 2 has focused on the interplay between the 
electronic susceptibility and the electron-phonon interac¬ 
tion as possible mechanisms for the CDW formation®®’®®. 
The large temperature dependence indicates that strong 
anharmonic effects survive over a pressure range that is 
way larger than the experimentally reported CDW sta- 
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FIG. 2. Phonon dispersion at 3.5 GPa (a) and 16 GPa 
(b). Closed (resp. open) dots represent the experi¬ 
mentally determined longitudinal acoustic (longitu¬ 
dinal optic) soft phonon dispersions, and lines repre¬ 
sent the results of phonon calculations. Close to the 
r point, the phonon dispersion shows clearly an op¬ 
tical branch. Closer to the M point, the energy of 
the phonon observed correspond to the acoustic like 
branch calculated. Colored thick and thin lines rep¬ 
resent the SSCHA results for the longitudinal acous¬ 
tic and optical anharmonic phonons, respectively, at 
several temperatures. The thin black lines represent 
other phonon modes present in this energy range cal¬ 
culated within the SSCHA at 20 K, which are very 
harmonic and barely depend on temperature (see SI). 


bility range (0 < P < Pcdw) and might therefore play 
an important role that was up to now largely disregarded. 

In Fig. 3, we plot the square of the soft phonon fre¬ 
quency measured at T = 20 K at qcow as a function 
of pressure. The linear dependence indicates that the 
mode hardens as ^/P — Pcdw- We also plot the square 
of the phonon frequency calculated within the harmonic 
approximation, which is negative - i.e. indicating that 
the system is unstable against CDW formation - up to 
14 GPa >> Pcdw- This is reminiscent of the situation 
encountered at ambient pressure in NbS 2 , where we have 
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FIG. 3. Pressure dependence of the square of the mea¬ 
sured energy of the longitudinal acoustic soft phonon 
at the CDW ordering wavevector and at 20 K (red 
dots) compared to the square of the energy calculated 
in the harmonic approximation (black diamonds). 
Lines are linear fits. The red region indicates the 
range of pressure where a CDW is observed exper¬ 
imentally at 20 K. The orange region indicates the 
range of pressure where harmonic calculations pre¬ 
dict a CDW ground state. 


observed a large phonon softening upon cooling, insuffi¬ 
cient however to induce the CDW formation despite the 
predictions of the harmonic ab-initio calculation^^. 

To quantify the effect of anharmonicity and check 
whether it is responsible for the destabilization of the 
CDW between Pcdw and 14 GPa, we have carried out 
a series of calculations for the pressures and tempera¬ 
tures experimentally investigated, within the SSCHA ap- 
proach^^’^^ (see SI) which allows us to access directly the 
anharmonic free energy of the system, with full inclusion 
of the anharmonic potential terms. 

The results are superimposed to our experimental dis¬ 
persions in Fig. 2. Not only is the instability suppressed as 
experimentally observed, but the agreement between the 
measured and calculated dispersions over the entire pres¬ 
sure range is remarkable compared to harmonic ab-initio 
calculation. A small overestimation of phonon calcula¬ 
tion is nevertheless observed in the range around qcDW- 

A key result here is that, even at low temperature, 
anharmonic terms play a key role in the destruction of 
the harmonically stable long-range CDW order. This is 
surprising at first glance since anharmonic effects usu¬ 
ally manifest themselves at finite temperature, through 
the effect of many-body interactions between thermally 
excited phonons. These account for the lattice thermal 
expansion and for the decay of the phonons lifetime at 
high temperatures, but are usually strongly suppressed at 
low temperatures. On the other hand, given that atoms 



FIG. 4. (a) Calculated phonon dispersion at 7 GPa 
and 0 K with the SSCHA. The size of the bars is pro¬ 
portional to the electron-phonon contribution to the 
phonon linewidth (the bar is twice the FWHM). (b) 
The Eliashberg function o?F(ijj) and the integrated 
electron-phonon coupling A(ti;).(c) The phonon den¬ 
sity of states decomposed into the different atoms. 
In all figures the grey background denotes the inte¬ 
gration area used to infer that the soft acoustic lon¬ 
gitudinal mode contributes 17% of the total electron- 
phonon coupling. 


vibrate even at zero temperature due to the zero-point 
fluctuations, their motion can in principle be affected by 
the anharmonic potential. This appears to be the case 
in this system where, below 14 GPa, the harmonic part 
of the potential that drives the lattice instability is over¬ 
whelmed by the anharmonic part. The balance between 
these two parts of the potential tells whether there is 
a lattice instability (driven by the CDW) or not. This 
bears stark analogy with the destruction of long range 
magnetic order from spin-spin quantum fluctuations in 
low dimensional quantum magnets^^. 

Finally we note that these many-body effects are rel¬ 
evant not only at q = qcDWj but extend up to the M 
point, the boundary of the Brillouin zone at q = a*/2, 
where two low-energy modes, the longitudinal acoustic 
mode that drives the CDW transition and an optical 
mode, are substantially hardened with respect to the re¬ 
sult of the harmonic calculation. 

The importance of the electron-phonon interaction in 
the CDW formation has been emphasized^^’^®, but its 
role in superconductivity remains debated, as is, more 
generally, the interplay between these two orders. Indeed, 
the absence of reliable ab initio phonon calculations in the 
high symmetry phase, has prevented a calculation of the 
superconducting critical temperature. Having achieved 
a detailed understanding of the low energy part of the 
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P (GPa) 

A 

ujiog (meV) Tc calc. (K) Tc 

exp. (K) 

3.5 

1.28 

12.6 

11.3 

7.8 

7 

1.16 

13.7 

10.5 

8.2 

16 

0.91 

16.4 

7.8 

7.8 


TABLE I. Calculated A, logarithmic frequency average oJiog, 
and Tc values for 2H-NbSe2 at 3.5, 7 and 16 GPa. The McMil¬ 
lan equation with n* = 0.16 is used to calculate Tc. Experi¬ 
mental data are taken from Ref.^°. 

experimental phonon spectra and having shown how an- 
harmonicity can suppress the long range ordered CDW 
phase, we can now address the superconducting prop¬ 
erties of NbSe 2 by calculating the Eliashberg function 
and the integrated electron-phonon coupling (see SI) as 
a function of pressure. As of now, this calculation can 
only be done in the high pressure phase, where the num¬ 
ber of atoms per unit cell is low enough. We show in Fig. 
4 the results of the calculation at 7 GPa (data for other 
pressures are shown in the SI). We find that the aver¬ 
age electron-phonon coupling is as large as A = 1.16 (see 
Tab. I for other pressures), demonstrating that NbSe 2 is 
a strong coupling superconductor. 

Previous works^^’^^’^®’^^’^^ identified the large 
electron-phonon coupling of the soft acoustic mode 
at qcDW as the mechanism responsible for the CDW 
formation. For this reason, it is natural to expect that 
this soft mode also plays a key role in superconductivity. 
However, its contribution to the average electron-phonon 
coupling is less than 17% (as it can be inferred from the 
integrated A(a;) in Fig. 4), meaning that it has only a 
marginal role in superconductivity. Indeed, the strong 
coupling of this mode is very localized around qcDWj 
and its contribution to A averages out when integrating 
over the Brillouin zone. The modes sustaining super¬ 
conductivity are given by the two main contributions to 
the Eliashberg function: a broad peak centred around 
« 11.5 meV and a second one in the 22 — 26 meV region. 
These two features contribute to « 70% of the total 
electron-phonon coupling. Further insights are obtained 
by decomposing the phonon density of states (DOS) 
in atomic vibrations along different directions. As the 
DOS and a^F{uj) are very similar, the decomposition 
will also apply to the Eliashberg function. The low 
energy feature at 11.5 meV is mainly due to Nb and 
Se vibrations parallel to the Nb plane, corresponding 


to the flat phonon band along PM suffering from a 
large anharmonic correction, while the features in the 
22 — 26 meV range are mostly attributed to out-of-plane 
and in-plane Se displacements with a non negligible Nb 
component. In both modes sustaining superconductivity 
there is a substantial Se component, upturning the 
conventional wisdom that in superconducting dichalco- 
genides the chalcogene is accessory while the transition 
metal plays a key role^. Our work demonstrates that 
the superconducting properties strongly depend on 
both the transition metal and the chalcogene, offering 
a natural explanation for the different relation between 
CDW and superconductivity encountered in different 
dichalcogenides^^^^"^. Moreover, as optical phonon modes 
contribute strongly to superconductivity in 2H-NbSe2, 
while the CDW is determined by the softening of the 
longitudinal acoustic mode, the superconducting Tc is 
insensitive to the presence of the CDW, explaining the 
phase diagram in Fig. 1. 

Our work points toward the key role of anharmonicity 
in destroying the CDW order not only through the usual 
thermal phonon-fluctuations, but also through quantum 
phonon-fluctuations, and on a very wide pressure range 
above the quantum phase transition. Our results, in this 
respect, are fundamental and relevant for a very large set 
of materials exhibiting a CDW instability - or more gen¬ 
erally a second order electronic quantum phase transition 
coupled to the lattice - such as transition metal dichalco- 
genides, cuprates, Bechgaard salts and transition metal 
bronzes. 
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